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We study the percolation properties of the growing clusters model. In this model, a number of 
seeds placed on random locations on a lattice are allowed to grow with a constant velocity to form 
clusters. When two or more clusters eventually touch each other they immediately stop their growth. 
The model exhibits a discontinuous transition for very low values of the seed concentration p and 
a second, non-trivial continuous phase transition for intermediate p values. Here we study in detail 
this continuous transition that separates a phase of finite clusters from a phase characterized by the 
presence of a giant component. Using finite size scaling and large scale Monte Carlo simulations we 
determine the value of the percolation threshold where the giant component first appears, and the 
critical exponents that characterize the transition. We find that the transition belongs to a different 
universality class from the standard percolation transition. 
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I. INTRODUCTION 



Percolation represents a paradigmatic model of a geo- 
metric phase transition [lHg ■ It has been widely studied 
and has numerous applications in many areas of physics 
[8l-fl7l|. Its importance and practical applications are de- 
scribed in detail elsewhere, see for example Q. Here, 
we will present for the sake of clarity and completeness, 
some necessary definitions of important quantities that 
commonly appear in the literature. In the classical site 
percolation model, the sites of a square lattice are ran- 
domly occupied with particles with probability p, or re- 
main empty with probability 1 —p. Neighboring occupied 
sites are considered to belong to the same cluster and per- 
colation theory simply deals with the number and prop- 
erties of these clusters. When the occupation probability 
p is small, the occupied sites are either isolated or form 
very small clusters. On the other hand, for large p there 
are a lot of occupied sites that have formed one large 
cluster and it is possible to find several paths of occupied 
sites which a walker can use to move from one side of the 
lattice to the other. In this latter case, it is said that a 
giant component of connected sites exists in the lattice. 
This component does not appear in a gradual "linear" 
way with increasing p. It appears suddenly at a critical 
occupation probability p c . Below p c there are only small 
clusters and even if we increase the lattice size consider- 
ably, these clusters remain small, i.e. the largest cluster 
does not depend on the system size. Above p c , suddenly, 
small clusters join together to form a large cluster whose 
size scales with system size. Hence, the term giant com- 
ponent or infinite cluster which is very common in the 
literature 

In percolation, p plays the same role as the tempera- 
ture in thermal phase transitions, i.e. that of the control 
parameter, while the order parameter is the probability 
Poo that a site belongs to the infinite cluster. For p > p c , 



Poo increases with p by a power law 

Poo ~ (P - Pcf 



(1) 



Other important quantities are the correlation length £ 
which is defined as the mean distance between two sites 
on the same finite cluster and the mean number of sites 
S of a finite cluster. When p approaches p c , £ increases 
as 



£ ~ (p-Pc)~ 



(2) 



The mean number of sites S of a finite cluster also di- 
verges at p c 



(3) 



The critical exponents /3, v and 7 describe the critical 
behavior associated with the percolation transition and 
are universal. They do not depend on the structure of 
the lattice (e.g., square or triangular) or on the type of 
percolation (site, bond or even continuum) [2j. 

In this paper we study numerically t he p ercolation 
properties of the growing clusters model [18j which we 
describe in Sec. [TTJ We focus on the intermediate con- 
centration regime and find that the model exhibits a non- 
trivial percolation transition which belongs in a different 
universality class from standard percolation. We deter- 
mine quite accurately the position of p c , and the values 
of the critical exponents associated with this transition. 



II. THE GROWING CLUSTERS MODEL 

The growing clusters model is presented in detail in 
[l8j |. Here, we provide a brief description of it. A square 
lattice of size L is randomly populated with "seeds" with 
probability p. These seeds are allowed to grow to clusters 
isotropically and stop when they touch another growing 
aggregate, see fig. [1] for a schematic of the system evo- 
lution. At every time step, i.e. one Monte Carlo Step 
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FIG. 1: Top Left: At time t = Monte Carlo Steps (MCS) 
6 seeds are randomly placed on a 40 x 40 lattice. Top Right: 
After t = 3 MCS there are 6 evolving clusters. Bottom Left: 
After t — 6 MCS there are 4 evolving clusters, as 2 clusters 
have touched each other and stopped growing. Bottom Right: 
Final state of the system. There are no evolving clusters and 
3 stable clusters have been formed. 

(MCS) all seeds are investigated once for the possibility 
to grow in size instantaneously in all neighboring sites. 
Investigation sequence is random in order. Each seed is 
allowed to grow its periphery by one layer (increase the 
radius by one) provided that there is no overlapping with 
other growing seeds. Thus, each seed becomes an evolv- 
ing cluster. As soon as two, or more, clusters touch each 
other, the growth of all of the adjoined clusters stop. 
Neighboring sites are considered to belong in the same 
stable cluster. To study the model we perform large scale 
Monte Carlo simulations as described in [l8[ and we mon- 
itor the properties of the clusters that are formed. 



III. FINITE SIZE SCALING 

Equations [TJ3] are valid for infinite systems close to the 
critical threshold. In practice, however, it is not possible 
to use them to calculate the critical exponents with con- 
siderable accuracy due to finite size effects. Thus, one 
has to resort to finite size scaling Q techniques. Due to 
the finite size of the lattices that can be simulated the or- 
der parameter is expected to depend on the system size. 
Assuming that we are close to the critical threshold so 
that the correlation length £ is comparable to the system 
size L it can be shown that the probability P ma x that 
a site belongs to the largest cluster follows the scaling 



relation: 

Pm m =L-^"F[L 1 ^(p-p e )] (4) 

where we have deliberately used the notation P ma x in- 
stead of Poo in order to emphasize the finiteness of the 
systems. Here F is a suitable scaling function. Simi- 
larly, any other quantity varying as p — p c \ x is expected 
to scale similar to eq. U with /3 replaced by the appro- 
priate exponent x and with a different scaling function 
F. At p = p c , these quantities are expected to scale as 
power law since the scaling function reduces to a sim- 
ple proportionality constant. This result gives a way to 
determine the critical exponents. One important char- 
acteristic of the standard percolation transition is that 
exactly at the critical point, the largest cluster has a 
fractal geometry meaning that the mass of the largest 
cluster Si scales with the system size as Si ~ L d s , where 
the fractal dimension df is known to be equal to 91/48. 
Since Si = P max L d , where d is the dimensionality of the 
space (d = 2 in the lattice case) one can easily derive a 
scaling law relating d/,/3 and u, namely dt = d — ft fv. 

IV. RESULTS AND DISCUSSION 

To obtain an indication of the system critical behav- 
ior, we start by monitoring the size of the largest cluster 
S-\ as a function of the initial "seed" probability p, see 
[18| and fig. 2 therein. There, it is evident that the 
system exhibits two phase transitions: One very sharp 
transition at p = and a second, smoother, transition 
around p ~ 0.5. The first transition is discussed in de- 
tail at [18[ and characterized by the fact that the size of 
the largest cluster Si, which is the order parameter of 
the system, has a discontinuity for p = and, thus, the 
system exhibits a sharp, albeit artificial, first order phase 
transition. 

The second phase transition turns out to be rather in- 
teresting. It is reminiscent of the classical percolation 
transition and it is important to clarify the extent of this 
similarity. Thus, we simulate systems of several different 
sizes and for several different initial seed concentrations. 
After allowing the growth process to complete and the 
systems to reach their final states, we study the proba- 
bility P m ax that a randomly chosen site belongs to the 
largest cluster. When there are only few initial seeds, af- 
ter the evolution of the system is completed, the clusters 
formed are small and isolated. However, we expect that 
with increasing p when a large number of seeds is intro- 
duced, the growth process will lead to the formation of a 
large spanning cluster. This is the classic behavior seen 
in a system which undergoes a continuous phase tran- 
sition. In such cases, we can use finite size scaling to 
determine the position of the critical concentration, p c , 
and the critical exponents. Initially, we are interested in 
the geometry of the largest cluster at criticality and its 
fractal dimension df. In fig. [5] we plot the size of the 
largest cluster Si — P m axL 2 as a function of the lattice 
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FIG. 2: Size of the largest cluster, Si = P mol L 2 , as a function 
of the lattice size L for three different initial concentrations, 
p = 0.48, 0.496 and 0.52(black squares, red dots and green 
triangles respectively). Points are Monte Carlo Simulation 
results and the straight line is a power law fit with slope 1.79. 



FIG. 3: Snapshot of the final state of a system with L = 100 
with initial seed concentration p = 0.496 which is equal to the 
critical concentration (p — p c ). Different colors correspond to 
different clusters. The final coverage is p/ = 0.538. The 
largest cluster is shown in black color. 



size L for three different initial concentrations, namely 
p = 0.48, 0.496 and 0.52 (black squares, red dots and 
green triangles respectively). For p = 0.48 (squares) we 
observe a downward bending curve which clearly indi- 
cates that p = 0.48 is below p c as Si does not scale with 
the system size for large L. For p = 0.496 (circles) we 
observe a straight line implying a power law scaling as ex- 
pected for p — p c . This value for the critical point agrees 
reasonably well with our preliminary estimation for the 
critical point from the position of the maximum of the 
second largest cluster [lg]- Forp = 0.52 (triangles) we ob- 
serve a curve which is bending slightly upwards for large 
L, indicating that we are above p Cl although admittedly 
for the sizes of lattices simulated here is quite difficult to 
observe the bending for p — 0.52. We will, however, use 
another method below to confirm that p c — 0.496 for the 
growing clusters model. We calculate the fractal dimen- 
sion, df, from the slope of the straight line at p c and find 
d f ~ 1.79. 

Figure [3] shows a snapshot of the clusters that have 
been formed in a random realization of a growing clus- 
ters system with L = 100 at p c = 0.496 once the system 
has reached its steady state. Periodic boundary condi- 
tions have been used for the cluster labeling. The largest 
cluster is shown in black color. 

In order to determine more accurately the percola- 
tion threshold and the critical exponent j3/v ratio in 
the same time, we use eq. |U We plot in fig. @] 
vs p for seven different system sizes, namely 
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FIG. 4: PmaxL vs p for seven different system sizes, 
namely L = 100, 200, 400, 600, 800, 1000, 1200. The curves 
cross at p c = 0.496 in agreement with the scaling relation 
eq. H 
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FIG. 5: Mean mass S of the finite clusters as a function of L 
for p = p c . Points are Monte Carlo Simulation results. The 
straight line is a best fit to the simulation data and has a 
slope of 1.63 



FIG. 6: PmaxL - 206 vs (p-p c )L ' 85 for seven different system 
sizes L = 100, 200, 400, 600, 800, 1000, 1200. Points are simu- 
lation data. The data collapse to one single curve. Thus we 
estimate 1/v = 0.85. 



L = 100,200,400,600,800,1000,1200, and we vary /3/u 
until all curves cross at one single point only. This is done 
for /3/f = 0.206 and for p = p c ~ 0.496. This result is in 
excellent agreement with our previous estimation for the 
df and the scaling relation df = d — @/v. It also allows 
to determine the exact location of the critical point with 
more accuracy than the previous method. 

Fig. [5] shows a plot of the mean mass S of all the 
finite clusters as a function of L for p = p c . At criticality, 
this quantity is expected to scale as S ~ L 1 ^ (see section 
III). From the slope of the straight line we estimate 7/1/ = 
1.63. This, as well as the previous result, are in very good 
agreement with the scaling relation dv = 7 4- 2/3. 

Finally, we calculate the critical exponent v. In fig. 
[5] we plot P max L - 206 vs (p — p c )L x / u for seven different 
system sizes L = 100,200,400,600,800,1000,1200, and 
we vary the exponent l/v until all data collapse on one 
single curve. The data collapse, in agreement with eqf?] 
enables us to determine the critical exponent v with con- 
siderable accuracy. We find that 1/u = 0.85. 

We can determine other critical exponents from the 
scaling relations that are known to connect them. Be- 
low, for completeness, we present a table with the critical 
exponents of the growing clusters model in comparison 
to those of classical percolation. The difference in the 
critical exponents shows that the two models belong in 
a different universality class. Moreover, our calculated 
value of v implies that the growing clusters model is also 
in a different universality class from the "invasion perco- 



lation" [TT[ model. 



Exponent 


Class. Percolation 


Growing Clusters 


P 


0.138 


0.24 


7 


2.38 


1.91 


V 


1.33 


1.17 


* = 1/(0 + 7) 


0.395 


0.46 


t — 1 + avd 


2.05 


2.08 


df 


1.89 


1.79 



TABLE I: Comparison of the critical exponents between the 
classical percolation transition and the transition of the grow- 
ing clusters model. The exponents r and a are associated with 
the cluster size distribution [l[ 



V. CONCLUSIONS 

We have studied the growing clusters model and found 
that it exhibits two phase transitions with increasing seed 
concentration. A first order transition at p = and a 
continuous transition at p c = 0.496, separating a phase 
of isolated clusters for p < p c from a phase where a giant 
component is present for p > p c . Using finite size scaling 
we have calculated the position of the phase transition 
and the critical exponents with considerable accuracy to 
establish that this transition belongs to a different uni- 
versality class from the standard percolation transition. 
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